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^ . Abstract 

The joint spectral radius of a finite set of real d x d matrices is defined to be the 
maximum possible exponential rate of growth of products of matrices drawn from 
that set. In previous work with K. G. Hare and J. Theys we showed that for a certain 
one-parameter family of pairs of matrices, this maximum possible rate of growth is 
\Q I attained along Sturmian sequences with a certain characteristic ratio which depends 

O ■ continuously upon the parameter. In this note we answer some open questions from 

that paper by showing that the dependence of the ratio function upon the parameter 
takes the form of a Devil's staircase. We show in particular that this Devil's staircase 
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Q ■ attains every rational value strictly between and 1 on some interval, and attains 

irrational values only in a set of Hausdorff dimension zero. This result generalises 
to include certain one-parameter families considered by other authors. We also give 
explicit formulas for the preimages of both rational and irrational numbers under the 

S^ . ratio function, thereby establishing a large family of pairs of matrices for which the 

joint spectral radius may be calculated exactly. 
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1 Introduction 

The spectral radius of a. d x d real matrix A, which we denote by p(y4), is defined to be 
the maximum of the moduli of the eigenvalues of A. If || ■ || is any norm on M'^, then 
the spectral radius satisfies the well-known identity p(A) = lim„_>oo ||A"||^/". Given a 
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bounded set A of real d x d matrices, we by analogy define the joint spectral radius of A 
to be the quantity 

^(A) := lim max |||Ai„ ■ ■ ■ A^J] - : A. G a| . 

n— ^co L ■'J 

It is not difficult to establish that this limit exists (essentially as a consequence of subaddi- 
tivity) and that its value is independent of the choice of norm || ■ ||. The joint spectral radius 
was introduced by G.-C. Rota and G. Strang in 1960 (see ll29l . later reprinted in ll28l ) and 
is the subject of ongoing research interest, which has dealt with its applications, its com- 
putation and approximation, and its intrinsic properties as a mathematical function. For a 
broad range of references on this topic we direct the reader to ll5l [T3l[T6ll26ll . 

It is not difficult to show that the joint spectral radius admits the alternative formulation 

j_ 
g{A) = sup limsup ll^n ■ ■ ■ y4i|| " , 

and that when A is compact there exists a sequence (Ai) of elements of A such that 
||y4„ ■ ■ ■ Ai||^/" — )■ f)(A). (A proof of this statement may be found in [fT6ll .') In this paper 
we are concerned with the following general question: given a finite set of matrices A 
and a sequence (Ai) in A such that \\An ■ ■ ■ Ai||^''" — )■ f?(A), what can we say about the 
structure of the sequence (A)? 

A question of particular interest is that of when there exist periodic sequences of ma- 
trices which achieve this maximal rate of growth. In ETI . J. Lagarias and Y. Wang asked 
whether every finite set A of rf x c? real matrices has the property that H^^ ■ ■ ■ Ai ||^/" — )• 
^(A) for some periodic sequence of elements of A, or, equivalently, whether every A has 
the property that g(A) = p{Ak ■ ■ ■ AiY^'' for some finite sequence Ai, . . . , A^ E A. We 
shall say that A has the finiteness property if such a periodic sequence exists. The existence 
of pairs of 2 x 2 matrices which do not satisfy the finiteness property was subsequently 
established by T. Bousch and J. Mairesse JH, with additional proofs being given later by 
V. Blondel, J. Theys and A. Vladimirov Q and V. Kozyakin [19]. The finiteness property 
continues to be the subject of research investigation: some sufficient conditions for the 
finiteness property have been given in ^dlHllTTl, and in a recent preprint N. Guglielmi 
and V. Protasov have given an algorithm for the rigorous verification of the finiteness 
property for real matrices lfT2l . 

In [fT3l , together with K. G. Hare and J. Theys the present authors investigated the 
finiteness property for pairs of matrices of the form Aq, := < Aq , A" L where 



and a G [0, 1]. It was shown in particular that if (xi) E {0, 1}^ is a sequence such that 
\\Ai^' ■ ■ ■ Axi \\'^^"' — )■ g{/Ka), then the proportion of terms of (xj) which are equal to 1 
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is well-defined and equal to a value x(a) E [0, 1] which depends only on a. We further 
showed that r is a continuous function of a, and gave an explicit expression for a value a^ 
such that r(Q;*) ^ Q, providing a completely explicit example of a pair of matrices which 
does not have the finiteness property (see formula (18.21 ) below). 

In this paper we undertake a detailed study of the behaviour of the function r for a 
belonging to the larger domain [0, oo). We extend the results described above in several 
directions. Firstly we give an explicit formula for t~^(7) when 7 G (0, 1) fl Q, and prove 
that this preimage is always an interval with nonempty interior. This allows us to construct 
an infinite family of examples of pairs of 2 x 2 matrices where the joint spectral radius 
may be computed exactly. Since the problem of devising algorithms for the computation 
of the joint spectral radius is ongoing (for some recent contributions see [fni5l [T2ll20ll26l ) 
these examples are potentially of value for the testing of new algorithms. 

Secondly, we show that the function r takes the form of a Devil's staircase, as was 
conjectured in [[T3l [30l . The methods which we use to obtain these first two results are 
significantly more general than those used in [fT3ll . and can also be applied to the families 
of pairs of matrices studied by other authors in [|4l [T9l|. We show in particular that r takes 
rational values only in the complement of a set of Hausdorff dimension zero. This result 
was previously noted in a special case in dH, though no proof was given. 

Finally, we give an explicit formula for t^^(7) when 7 G (0, 1) \ Q, and in the special 
case of the matrices given by (11.11 ) we provide some inequalities for use in the rigorous 
computation of its value. We thus show how to construct an uncountable family of explicit 
examples for which the finiteness property is not satisfied. As with our explicit description 
of pairs of matrices which satisfy the finiteness property, we anticipate that these examples 
may be of value in future in the analysis of algorithms for computing the joint spectral 
radius. 

2 Notation and statement of results 

Throughout this paper we will consider pairs of real 2x2 matrices which we denote by Aq, 
Ai. To describe the structure of sequences of these matrices we use the space of symbolic 
sequences E2 := {0, 1}^. We refer to the elements of S2 as infinite words. We equip S2 
with the infinite product topology, with respect to which it is compact and metrisable. On 
some occasions it will be useful to employ a metric on S2: to this end, given sequences 
{xi),{yi) G S2 we define 

ci[(x,),(y,)]:=2---i---'=^'>, 

where the expression 2~°° is interpreted to mean 0. This defines an ultrametric on E2 
which generates the infinite product topology. We also define the shift transformation 
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T: S2 — 7- S2 by T[(xj)] := (xj+i), which is a continuous surjection. If a pair of matrices 
B := {Bq, Bi} is understood, then following the terminology of |fT3ll24ll we shall say that 
a sequence x = (xj) G S2 is weakly extremal for B if \\Bx^ ■ ■ -E^JI^/" — > g{B) in the 
limit as 72 —i- 00. 

In addition to considering infinite sequences in {0, 1} we shall also find it useful to 
consider finite sequences, which we refer to as finite words. If m = (ui)"=i is a finite word 
we call n the length of u and define \u\ := n.To simplify certain statements we allow the 
word of length zero, which we refer to as the empty word. 

With the pair of matrices A = {Aq, Ai\ fixed, we define A^^' := Aq and A]^' := aAi 
for all real numbers a > 0, and let Aq, := {Aq , Ai"'}. We shall denote the quantity 
g{/\a) simply by Q{a). The function g: [0, 00) — > M is continuous (see for example 
[fT4ll ). For every x G S2, n > 1 and a > we define Aa{x,n) := Ax^ ■ ■ ■ Ax^ and 
A{x., n) := Ai{x., n) = Ax„ ■ ■ ■ Ax^. If m is a finite word of length m > 1, we similarly 
define A^{u) = Ai"i ■ ■ ■ ^1°^ and A{u) = Ai{u). 

In this paper we are concerned specifically with pairs of matrices such that the max- 
imum growth rate of partial products occurs along Sturmian sequences. A large range 
of definitions of Sturmian sequence exist in the literature, see for example [12511 and the 
surveys in IIT0ll22l . The definition which we give in this section is not the most straight- 
forward to state, but is the most suited to the proof methods which are used later in this 
article. In order to state this definition and describe its main consequences, we require 
some further terminology. 

Given a finite word u, let |m|i denote the number of entries of u which are equal to 1, 
and if w is not the empty word, define the slope ofu to be the quantity q{u) := |m|i/|m|. If 
u = (Mi)"=i and V = (fj)™i are finite words then we define the concatenation ofu with 
V, denoted by uv, to be the finite word co = (wj)"^™ such that cui = Ui for I < i < n and 
cui = Vi-n for n < i < n + m. We use the symbols and 1 to denote the words of unit 
length with entries and 1 respectively. For positive integers k we use the notation u'' to 
denote the successive concatenation of k copies of the word u, and we define m° to be the 
empty word. The word u'' will be referred to as the fcth power of u. Using these notational 
conventions it is clear that any finite word may be written in the form l'^'=0"'=-i ■ ■ ■l'^^ 
for some finite collection of non-negative integers Oj. Given a finite word u of nonzero 
length n, we use the symbol m°^ to denote the unique infinite word x = {xi)°Zi ^^ch that 
Xi+kn = Ui for all A; > and 1 < i < n. 

We say that the finite word m is a subword of the finite word v if v = aub for some 
(possibly empty) finite words a and b. If a is empty then we say that u prefixes v. We shall 
also say that a finite word u prefixes an infinite word x E T.2 if Ui = Xi for all i in the 
range 1 < i < \u\. A word u will be called balanced if for every pair of subwords vi, 
V2 of u with |fi| = |t>2| we have ||fi|i — |f2|i| < 1. Clearly u is balanced if and only 
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if every subword of u is balanced. We say that a; G S2 is balanced if every prefix of x 
is balanced. We say that two finite words u = (ui), v = (vi) are cyclically equivalent 
if they are equivalent by some cyclic permutation: that is, they share same length n and 
there exists an integer k such that Ui = Vi^k for 1 < i < n — k and Mj = Vi^k-n for 
n — k < i < n. It is not difficult to see that u and v are cyclically equivalent if and only 
if there exist (possibly empty) finite words a and b such that u = ah and v = ha. We say 
that u is cyclically balanced if it is balanced and all of its cyclic permutations are also 
balanced. One may show that a nonempty finite word u is cyclically balanced if and only 
if u°° is balanced (see e.g. [[131 Lemma 4.7]). 

An infinite word x G S2 will be called Sturmian if it is balanced and recurrent with 
respect to T. It follows that if m is a finite nonempty word, then m°° is Sturmian if and 
only if u is cyclically balanced. The key properties of Sturmian sequences are outlined by 
the following theorem, the proof of which may be found in [|22l[25l . 

Theorem 2.1. For each 7 G [0, 1] define a set X^ C S2 as follows: we have x G X^ if 
and only if there exists 5 G M such that either 

Xn = \j{n + 1) + (5J - \jn + 5\ 

or 

Xn = \l{n + 1) + 6] - \-fn + 6] , 

where [x\ = max{n G Z : n < x} and [x] = niin{r2 G Z : n > x}. Then an infinite 
word X G S2 is Sturmian if and only if x E \J^^\q 1] X^- The sets X^ have the following 
properties: 

(i) Each X^ is compact and satisfies TX^ = X^. 

(ii) The restriction ofT to X^ is uniquely ergodic, i.e., T has a unique invariant mea- 
sure. 

(Hi) Ifx G Xy then n^^^{l < i < n: x j = 1 } — )■ 7 a^ n — )■ cxd. 

(iv) Ifj= p/q in least terms then the cardinality of X^ is equal toq.If'j is irrational 
then X^ is uncountable. 

Whilst our primary objective is to continue the study of the pair of matrices defined by 
(11.11) which were examined in [[3l[T3l[30l, the methods which we use are general enough 
to encompass a larger family. The following definition describes the minimum properties 
necessary for our arguments to apply: 

Definition 2.2. Let A = {Aq, Ai} be a pair of 2 x 2 real matrices. We shall say that A 
satisfies the technical hypotheses if the following properties hold: 
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(i) The matrices Aq and Ai are non-negative, invertible, have positive trace, and do not 
have a common invariant subspace. 

(ii) If M is a finite word which is not of the form 1" or 0" then all of the entries of the 
matrix A{u) are positive. 

We shall further say that A satisfies the Sturmian hypothesis if there exists a function 
r : [0, oo) -T- [0, 1] such that the following properties hold: 

(iii) For each a > 0, every x G X^(a) is weakly extremal for A^. 

(iv) For each a > 0, if x e S2 is weakly extremal for Aq, then r7,~^#{l < i < n: Xi = 
1} ^r(a). 

(v) If n is a finite word which is not cyclically balanced then p{Aa{u)) < ^(ct)'"'. 

The function t will be called the 1-ratio function of the pair A. 

Note that as a consequence of (iv), if r exists then it is unique. By Il24l Theorem 2.3] 
and the minimality of the invariant sets X^, the hypothesis (iii) is in fact equivalent to 
the hypothesis that X^(^a) contains at least one extremal infinite word. Some conditions 
equivalent to (iv) have been used by other authors: a description of these conditions and a 
proof of their equivalence are given in fl24l §6]. 

A range of examples of pairs A which satisfy the Sturmian hypothesis are known. 
In lfT3l . the authors together with K. G. Hare and J. Theys proved that the family of 
matrices given by (11.11) satisfies parts (iii)-(v) of the Sturmian hypothesis for a restricted 
to the interval [0, 1]. If we extend the definition of r to the interval [0, 00) by defining 
r(a) = l — x{l/a) foreacha G (1, 00), then by taking advantage of the relation Aq = Af 
it is not difficult to show that the Sturmian hypothesis in full for the family A^. The 
essential points of this argument are contained in Lemma [3711 below. 

In the earlier work ||4l, T. Bousch and J. Mairesse also proved that the Sturmian hy- 
pothesis holds for the matrices 

Ao:=[ ^. J, M:=[^ ^.,,^J, (2.1) 

subject to the inequalities K,hQ,hi > and Hq + hi < 2. Clearly the examples given by 
(|l.ll) and (12.11) also satisfy the technical hypotheses. In a series of papers, V. S. Kozyakin 
has shown that the Sturmian hypothesis holds for pairs of triangular matrices having the 

form 

fa b\ , A 

^° ■= [0 1) ' ^^ ■= [c d 
where < a, c? < 1 < 6c; an overview of this work is given in [fT9l . Note that the 
examples considered by Kozyakin satisfy the technical hypotheses in the case b,c > 0, 
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and are simultaneously similar to a pair of matrices satisfying the technical hypotheses 
when 6, c < 0. 

In order to state the explicit formula for the intervals x^^{p/q) which forms part of 
our first theorem we require one last definition, namely that of a standard pair. The set 
of all standard pairs, which we denote by V, is defined to be the smallest nonempty set 
of ordered pairs of finite words which satisfies the following two properties: (0, 1) G V, 
and if (m, v) eV then {uv, v) E V and (m, vu) G V. We say that cj is a standard word if 
it is one half of a standard pair. Every standard word is balanced (see, e.g., [22, Proposi- 
tion 2.2.15]). If {u, v) is a standard pair then (u, vu^) and (wf", v) are also standard pairs 
for every n > 0, and it follows that every power of a standard word is a subword of some 
standard word. In particular every power of a standard word is balanced, and consequently 
every standard word is cyclically balanced. A highly detailed analysis of the properties of 
the set V may be found in Il22ll . 

The main result of the present paper is the following theorem: 

Theorem 2.3. Let Aq, Ax he a pair of! x 2 real matrices which satisfies both the technical 
hypotheses and the Sturmian hypothesis. Then: 

(i) The function r is continuous and monotone non-decreasing, and satisfies r(0) = 
and limQ,_>oo t(a) = 1. 

(ii) For every rational number 7 G (0, 1), the set t^^(7) is a closed interval with 
nonempty interior The interval t^^(O) fl (0, 00) is nonempty if and only if Aq is 
diagonalisable, and similarly t^^(l) is nonempty if and only if Ai is diagonalis- 
able. 

The intervals t^^{p/q) may be computed exactly by the following procedure. If Aq 
is diagonalisable, let Pq := lim„^ooP(^)""^o- Then 



'\0) 



0, 



L ' PiPoAi)\' 

Similarly ifAi is diagonalisable and Pi := lim„_^oo p{Ai)~"-A", then 

'p{PiAo) 



t-Hl) 



-00 



If p/q G (0, 1) in least terms, then there exists a standard pair {u,v) such that 
<,{uv) = p/q. Define \u\ := gi and \v\ := q2, let Bi := A{u), B2 := A{y) and 
A := B1B2, and let P := lini„_j.ooP(^)~"^" be the Perron projection associated 
to the positive matrix A. Then we have 

and g{a) = p{Aa{uv)Y^'^ for all a G t^^ip/q)- 
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(in) The Hausdorff dimension of the sett ^([0, 1) \ Q) is zero. In particular, r ^(7) is a 
singleton for any irrational 7 G (0, 1). 



Some examples of the explicit formulae generated by Theorem I2.3r ii) are given in 
Table |2.ll Note that a direct consequence of Theorem (23^11) is that when all of the entries 
of Aq and Ai are rational, the endpoints of x^^{p/q) are algebraic numbers of degree 
either 1 or 2. In the latter case both endpoints belong to the same quadratic field. 

As a direct consequence of (iii) we obtain the following result: 

Corollary 2.4. The function x is not Holder continuous. 

To illustrate the behaviour of the function r we reproduce a diagram from [13|: see 
Figure [T] below. 




Figure 1: This figure shows the graph of r for a restricted to the interval [0,5/4] for 
the family of matrices given by (11.11) . Using the explicit formula for t~^{l/n) given in 
Table 12. II one may show that in this case a~^x{a) — )■ 1/e in the limit as a — )■ 0. 



In the cases studied in ||4l |T3l |T9l| the continuity of the function r is established by 
using the particular characteristics of the matrices Aq, Ai in quite a strong fashion. In 
proving part (i) of Theorem 12.31 we observe that the continuity of the function r is in fact 
a corollary of its defining properties. T. Bousch and J. Mairesse have asserted in ||4l that 
part (iii) of Theorem 12.31 holds for the case of triangular matrices of the form (|2.1I) . but 
their proof remains unpublished. In |fT3l we proved for the case of matrices (11.11) that 
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1 -ratio 7 


Standard pair 


Interval r ^(7) 


1/2 


(0,1) 


4 5 
.5' 4. 


3/7 


(00101,01) 




(5+i§v^)7 (13+2^42)2 
(13+2v^)5 ' (| + ^v^)7 




2/5 


(001,01) 




(2+i^v'6)^ (5+2^6)2 
(5+2^6)3 ' (| + iiV6)5 




1/3 


(0,01) 




69-16v^ 1656-384^3 
72 ' 1331 




2/7 


(0001,001) 




(^ + %^^)' (9+4v^)3 
(9+4^5)4 ' (2+|iV5)7 




1/4 


(0,001) 




496-64^21 13671-1764^21 
441 ' 10000 




1/5 


(0,0001) 




10612-5261v^ 43466752-21549056v^2 
8192 ' 28629151 




1/6 


(0,00001) 




l + f v/5 ' (3+f V5)6 




1/7 


(0,000001) 




4+^15 ' i'i+^y 




1 

n+l 


(0,0"-il) 






N 71 + 1 






7 (l+t + W"'+4")" 


V 2 2ri2+8n 7 



Table 2.1: This table gives some examples of explicit formulae for the intervals ^"^(7) 
for the family of pairs of matrices given by (11.11) . Note that the endpoints of r^^(l/r2) are 
asymptotically equal to e/n + o{l/n) in the limit as 77, — )■ 00, a feature which may be 
observed in Figure [IJ 
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t~^(7) is a singleton if 7 is irrational and not Liouville. Theorem 12 .3^ 11) shows that this 
remains true for 7 irrational and Liouville. 

In the course of proving part (iii) of Theorem l2.3l we are able to establish the following 
result: if L C (0, 1) is a compact interval, then there exist constants K > 1 and 9 E (0,1) 
depending on L such that the interval t^^{p/q) has diameter less than K9'' for every 
p/q E L, where the fraction p/q is understood to be given in least terms. Heuristically, 
this result tells us not only that values of a for which r(a) is irrational are extremely 
scarce, but also that values for which x:{a) is a rational number with large denominator 
are still relatively scarce, at least when a lies within a given neighbourhood bounded away 
from zero and infinity. This result is given as Corollary 16 . 61 below. 

The second result of this paper is the following theorem which gives an infinite product 
formula for r^^(7). 

Theorem 2.5. Let Aq, A\ he a pair of matrices which satisfy the technical hypotheses 
and the Sturmian hypothesis, and let 7 E (0, 1) \ Q. Let (a„)^i E N^ be the sequence 
of continued fraction coefficients of'-y, and for each n > 1 let pn/qn be the corresponding 
convergent. Define a sequence of finite words (s„) inductively by setting s_i := 1, Sq := 0, 
Si := Sq'^^s^i and s„+i := s'^+^Sn-ifor every n > 1, and for each integer n > —1 
define Pn := p{A{Sn))- Then t^^(7) is the singleton set whose unique element is given by 

a, :=Um(^ = ^ TT ( ^" • (2-2) 

In general it is not clear whether the infinite product given here will always converge 
unconditionally, although we are able to prove this in special cases. For the family of 
matrices defined by (II. 1|) we are able to give a checkable criterion for a rigorous bound 
on the error in approximating a^ by partial products of the infinite product given above. 
The details of these estimates are given in §8. 

3 Convex analysis and continuity of the 1 -ratio 

In this section we give the proof of part (i) of Theorem 12. 3[ and introduce a concave 
function S: [0, 1] — )■ M which characterises the rate of growth of A along Sturmian tra- 
jectories. We begin with the following simple lemma: 

Lemma 3.1. Let A = {Aq, Ai} be apair of matrices which satisfies the Sturmian hypoth- 
esis, and let t: [0, oo) — )■ [0, 1] be the corresponding 1 -ratio function. Define a new pair 
of matrices A := {Aq, Ai} by Aq := Ai, Ai := Aq. Then A also satisfies the Sturmian 
hypothesis, and if x denotes the 1-ratio function of A, then x{a) = 1 —x{l/a) for all 
a E (0, oo). 
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Proof. Let us define Aq, := {Aq, aAi} for each a > similarly to the definition of A^,, 
and let ^{a) = g{Aa) for all a > 0. Define Aa{x, n) = AxJ ■ ■ ■ Axi for all x G S2 and 
n > 1, and for each x = (xj) G S2 define a new sequence x by Ti := 1 — Xi. We have 
X G X^ if and only if x G Xi_^. Finally, define'?: [0, 00) — )• [0, 1] hyx{a) := 1 — t(l/a) 
for all a G (0, 00), and 'c(O) := 0. Note that for all x G S2, ra > 1 and a G (0, 00) 
there holds the identity Aa{x,n) = a'^Ai/a{x,n). As a direct consequence we have 
g{a) = ag{l/a) for all a G (0, 00). 

We may now verify directly that A satisfies the Sturmian hypothesis with t" being its 1- 
ratio function. The case a = being trivial, let us fix a > 0. If x G X^a) , then x G X^(i/a) 
and therefore 



lim 



t/T./-!/ IX, /t I 



V" ._ M Ml/n 



lim q;"^i/q,(x, n) = 0^(1/0;) = ^(a) 
as required. If Hm^^oo ||^a(2;, n) || ^/" = ^(a) then by the same token we have 

U/c 



lim ||^i/„(x,n)f /" = ^(1/a) 



and therefore 

lim —4j^ {1 <Xj <n: x,- = 1} = lim ( 1 7^ {1 <Xj <n: x,- = 1} ) 

n->-oo n n->oo \^ ?2 / 

= l-r(l/a) =r(a) 

since Definition 12. 2r iv) applies to Ai/q,. Finally, if m = (Mj)i=i is a finite word which is 
not cyclically balanced, then the finite word u = (mj) -=1 defined by Mj := 1 — Mj is clearly 
also not cyclically balanced and hence 

p(X(m)) = «l"lp(A/a(M)) < al"lf?(l/«)l"l = Q{ap 
as required. The proof is complete. D 

The following general theorem was proved in [|24l : 



Theorem 3.2. Let Abe a metric space and let Aq, Ai : A — )■ M2(M) be continuous func- 
tions such that Ao(A) ^ Ai{\)for all A G A. Suppose that there exists a function r : A — )■ 
[0, 1] with the following property: for every x G S2 such that ||Aa.^(A) ■ ■ ■ A3;^(A)||^/" — )■ 
^({Ao(A), Ai(A)}) in the limit as n ^ 00, we have n~^{l < z < ra: Xj = 1} — > t(A). 
Then the function r is continuous. 

We may now directly deduce several parts of Theorem 12. Sf i). 

Lemma 3.3. Let A be as in Theorem \2.3\ Then the 1 -ratio function x: [0, 00) — )■ [0, 1] is 
continuous and satisfies r(0) = and limQ,_5.oo t(a;) = 1. 
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Proof. The continuity of r follows immediately from Theorem 13 .21 In the case a = it is 
obvious that ||^Q,(a;, n)||^/" -)■ p{Ao) > when x E Xq and ||^Q,(a;,n)||^/" -)■ for all 
other X, and it follows that c(0) = 0. In particular we have limQ,_^o ^(ct) = by continuity. 
Let A and r be as in Lemma 13.11 applying the preceding arguments to A it follows that 
liuia^Q't{a) = 0, and therefore 



lim r(a) = lim (1 — r(l/a;)) = Um (1 — r(a)) = 1. 



D 



The following proposition, which characterises t in terms of a concave function on the 
unit interval, forms the cornerstone of the proof of Theorem |23l In the special case where 
Aq and Ai are as defined by (11.11) . the results of Proposition l3.4l correspond approximately 
to those of [fT3l, Proposition 6.1]. 

Proposition 3.4. Let A := {Aq, Ai} be as in Theorem \2.3\ Then there exists a continuous 
concave function S: [0, 1] — t- M with the following properties: 

(i) For each 7 G [0, 1] we have ^ log ||^(x, n) \\ — )• 5(7) for every x G X^. 

(it) For each a E [0, 00) and 7 E [0, 1] we have e^^'^^a"' < g{a) with equality if and 
only if 'y = r(a). Consequently, for nonzero a we have r(«) = 7 if and only if 
— log a is a subgradient ofS at 7. 

(Hi) If u is a word of length kq with <j('u) = p/q, then {kq)~^ logp(^(-u)) < S{p/q), 
with equality if and only ifu is cyclically balanced. 

Proof. We will show first that there exists a function S: [0,1] — )• M such that properties 
(i)-(iii) hold, and show only at the end of the proof that this function is continuous and 
concave. We begin by constructing a function S which satisfies (i). Since Xq contains 
only the single point 0°° it is clear that (i) holds for 7 = with S'(O) := logp(y4o), and 
similarly for ^(l) := logp(y4i). Let us therefore consider 7 E (0, 1) and x E X^. Using 
Lemma [3J] we may choose a > such that r(a) = 7. For each n > 1 we have 

log||^(x,n)|| =\og\\Aa{x,n)\\ -#{1 < j < n: Xj = 1} -loga 

and it follows by Definition 12. 2r iii) together with Theorem 12. If iii) that 

lim — log ||^(a;, n)|| = log^(Q:) — 7loga. (3.1) 

n— 5>oo Tl 

Since the left hand side of this equation does not depend on the choice of a G t~^(7) and 
the right hand side does not depend on the choice of x G X^, we conclude that the identity 
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(I3.1|) holds for all such choices. In particular if we define S{'j) := log g{a) — 7 log a then 
(i) is satisfied. Now, if a G (0, 00) is given, then for any 7 G [0, 1] we have for all x G X^ 

5(7) + 7 log a = lim — (log ||^(a;,n)|| + # {1 < j < n: x,- = 1} ■ log a) 

n— >oo Ti 

= lim — log ||^Q,(x, n)|| <log^(a). 

It follows from parts (iii) and (iv) of Definition [2]2] that the above inequality is an equality 
if and only if 7 = x{a). This proves (ii) for all cases except those in which a = 0, which 
it is trivial to verify directly. 

Let us now prove (iii). Let u be as in the statement of the Proposition, and define 
X := u°^ G T.2. If u is cyclically balanced then x is Sturmian, and since q{u) = p/q we 
necessarily have x G Xp/g. Using (i) together with Gelfand's formula we may obtain 

S{p/q) = lim ——log\\A{x,nkq)\\ = lim — — log ||^(u)"|| = — log p{A{u)). 

TirCq Tl — ^00 TlrCq rCq 



n—^oo 



Now suppose that u is not cyclically balanced. Using Lemma 13.31 let us choose a > 
such that x{a) = p/q. Using Definition 12. 2r v) together with part (ii) above we obtain 

p{Aiu)) ■ a'^P = p{AM) < g{a)^'' = e'^^'^^P/^^a^P 

and therefore {kqY^ log p{A{u)) < S{p/q) as required to prove (iii). 

It remains to show that S is continuous and concave. As a consequence of (ii) we have 
for every 7 G (0, 1) 

5(7)= inf log f)(a) — 7 log a. 

ae(0,oo) 

The restriction of S to (0, 1) is thus an infimum over a set of affine functions of 7, and 
hence is concave. It follows from standard results in convex analysis that the restriction 
of S to (0, 1) is also continuous. 

Finally let us show that S is continuous on [0, 1], and hence is also concave on that 
interval. We will show that S is continuous at 0, the case of continuity at 1 being similar. 
Since S is concave, the limit of 5 at exists, so it suffices to show that there exists a 
single sequence (7„) of elements of (0, 1) which converges to zero and has the property 
that 5'(7„) converges to 5(0). To this end, let us choose a strictly increasing sequence of 
integers (rij) such that the sequence ||Aq^||^^Aq^ converges to some matrix P. For each 
n > it is not difficult to see that the word 0"1 is cyclically balanced, and therefore 
5'(l/(n+l)) = (n + l)"Mogp(A;jAi) using (iii). We thus have 

lim S{l/{n, + 1)) - S{0) = lim ^- logp{A^o'Ai) - logp (Aq) 

j->oo j->oo rij + 1 ^ 

= lim ^- (iogp(A;;Mi) -iog||A;;^||) 

J-5-00 rij + 1 " "^ 



lim^-logp(||A^^||-'A^Mi) 
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since p(mQ^ || Aq^Ai) converges to p{PAi) as j — > oo. The proof is complete. D 

The following result together with Lemma [33] completes the proof of Theorem 1231 1) . 

Corollary 3.5. The 1-ratio function x is non-decreasing. 

Proof. It is an elementary fact in convex analysis that if Ai and A 2 are subgradients of a 
concave function at 71 and 72 respectively, and 71 < 72, then Ai > A2. The result now 
follows by Proposition 13 .41 D 

4 Standard words 

In this section we exploit some well-known features of Sturmian words to obtain a pair of 
propositions dealing with the combinatorial structure of the sets X^jq. 

Let us define two maps from the set of standard pairs V to itself by r(M, v) := (m, uv) 
and A(n,t>) := {vu,v). Throughout the remainder of the paper we use the following 
notation for continued fractions: if ai , . . . , a„ are positive integers, then we use the symbol 
[ai, . . . , a„] to denote the finite continued fraction 



Pn 


i 


Qn ' 


1 
^ 1 




"i 1 

a2 + ... + — 
an 



where p„ and g„ are coprime. Given positive integers ai, . . . , a„ let us write pk/qk '■= 
[oi, . . . , Ofc] for all k in the range 1 < k < n, and define also po := 0, go := 1^ P-i '■= !» 
g_i := 0. Subject to these conventions, for each integer k in the range 1 < k < n the 
integers p^, q^ satisfy the recurrence relations pk = akPk-i +Pk~2 and g^ = ctfcQ'fc-i + q'fc-2- 
If {ai)°Zi i^ sequence of positive integers, we use the notation 7 = [ai, 02, . . .] to mean 

lim„^oo[ai, ••-,««] = 7- 

The following proposition will be applied in the proof that r^^([0, 1] \ Q) has zero 
Hausdorff dimension: 

Proposition 4.1. Let 7 = Pn/qn = [ai, • • • , On] cind Pn-i/qn-i = [cti, • • • , ctn-i] where 
n > 2 and a„ > 1. Then there exist x G X^ and an integer k > ^qn such that 

d{x,Ti"-^x) < 2-^ andq-\\ogp{A{T^x,qn-i)) = S{pn-i/qn-i). 

Proof. Let s_i = 1, sq = and Si := Sq^^^s_i. Define s^ inductively for 1 < /c < n by 
Sk '■= sl'liSk-2- For fc = 1 we have (sq, si) = r"i"^((0, 1)). An easy proof by induction 
shows that for odd k > 1, 

(sk-i, Sk) = (r«'= o A'^'^-i o . . . o A-^^ o ri-^)((o, 1)) 
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and for even k, 

{sk, Sk^i) = (A"'^ o ^'=-1 o . . . o A'^^ o ri-i)((0, 1)), 

so in particular each Sk is standard, and hence is cyclically balanced. Define Pk/qk '■= 
[ai, . . . ,ak] for 1 < k < n. We have |si|i = 1 = pi, \si\ = ai = qi, and \sk\i = 
otfc|sfc-i|i + |sa;_2|i and \sk\ = cbk\sk-i\ + \sk-2\ for 1 < A; < n. It follows by induction 
that \sk\i = Pk and |sfc| = g^ for 1 < A; < n. In particular we have? (s„) = Pn/qn = 1, and 
since s„ is cyclically balanced it follows that x := s'^ & X^. The formula s„ = s'^_^Sn-2 
implies that the infinite word x is prefixed by the finite word s°^_^ and the infinite word 
j^gn-ij, ig prefixed by the finite word s^l^\ In particular we have d{x,T^"-^x) < 2^^ 
where k = {an — l)qn-i- We have 

qn = an^n-i + g„-2 < («„ + l)g„-i < 3(a„ - l)g„_i, 

since a„ > 1, and thus k > |g„ as claimed. Finally we note that T^x is prefixed by the 
finite cyclically balanced word s„_i and so by Proposition [3]4tiii) we have 

q~\\ogp{A{T''x,qn-i)) = |s„_i|~Mogp(^(s„_i)) = S {pn-i / qn-i) ■ 

The proof is complete. D 

The following proposition will be used in the proof that t~^{p/q) is an interval: 

Proposition 4.2. Let p/q G (0, 1) with p and q coprime. Then there exists a standard 
pair {u,v) such that <,{uv) = p/q. If {u,v) is such a pair then \uv\i = p, \uv\ = q, 
\u\ ■ \v\i — ImIi ■ |f I = 1, the words {uv)"'u and {yu)^v are cyclically balanced for all 
n > 0, and the word v^v"^ is not cyclically balanced. 

Proof. Let us write p/q = [ai, . . . , a„] with a„ > 1. If n = 1 then define {u, v) := 
r""-2(0, 1). For odd n > 1, define 

{u, v) := (r--^ o A""-^ o . . . o A"^ o r"^-^)((0, 1)), 

and for even n define 

{u, v) := (A""-^ o r""-i o . . . o A"^ o r"^'i)((o, 1)). 

A proof by induction on n similar to that in the previous proposition shows that \uv\i = p 
and \uv\ = q, and hence there exists a standard pair (u, v) such that <j(Mf ) = p/q. 

For the rest of the proof we let {u, v) be such a standard pair. Since {u, v) is standard, 
it follows by definition that the pairs (A" o T){u, v) = ({uv)^u, uv) and (F" o A)(m, v) = 
{vu, (t>u)"f ) are standard pairs for all n > 0. In particular, the words {uv)"'u and {vu)'^v 
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are standard for all n > 0, and hence these words are cyclically balanced. It is easy to see 
that the set of all standard pairs (a, b) such that |a| ■ |6|i — |a|i • |6| = 1 contains the pair 
(0, 1) and is closed under the action of T and A, and it follows that every standard pair 
has this property. In particular we have ImI-IiiIi — |M|i-|t>| = las claimed. Furthermore, 
since {u,uv) is a standard pair we have |M|(|ut>|i) — |n|i(|Mf |) = 1 so that \uv\ is coprime 
to \uv\i, and so any standard pair {u,v) which satisfies ^(wf) = p/q necessarily has 
\uv\i = p, \uv\ = g as claimed. 

Finally, an easy inductive proof starting with the pair (0, 1) shows that for every stan- 
dard pair (a, b) there is a (possibly empty) finite word p such that ab = pOl and ba = plO, 
see [|22l p. 57]. Since {u,v) is a standard pair we have uv'^u = pOlplO for some fi- 
nite word p, and hence in particular u'^v'^ is cyclically equivalent to a word of the form 
OpOlpl. Since |lpl|i = 2 + |OpO|i the word OpOlpl is not balanced, and hence n^f^ is not 
cyclically balanced. D 

Remark 4.3. It is worth noting that if we put pi = |M|i,gi = \u\ and p2 = |^|i5Q'2 = |^|, 
then pi/qi and p2/q2 are the Farey parents of p/q: that is, they are the unique fractions 
such that pi/qi < p/q < P2/q2, < qi, q2 < q, qi + q2 = q and pi + P2 = P- In fact 
one may show that the pair (m, v) specified by Proposition 14.21 is unique, but this is not 
required for our argument. 

Example 4.4. Let p/q = 3/7; then u = 00101, u = 01. In particular we have u'^v'^ = 
00101001010101, which contains the subwords 0010100 and 1010101 and thus is not 
balanced (let alone cyclically balanced). 

5 Preimages of rational points 

In this section we give the proofs of the various clauses of part (ii) of Theorem l2.3[ Since 
S: [0, 1] — )■ M is concave it follows from elementary convex analysis that the left deriva- 
tive 5^(7) and the right derivative 5^(7) both exist and are finite for every 7 G (0, 1), the 
right derivative 5*^(0) at either exists or equals +00, and the left derivative ^'^(l) at 1 
either exists or equals —00. Furthermore, the set of all subderivatives of S* at 7 G (0,1) 
is precisely [>S'^(7), >S'^(7)], the set of subderivatives of S at is precisely (—00, 5'^(0)], 
and the set of subderivatives of 5 at 1 is precisely [5*^(1), +00). Note that the latter two 
intervals are empty if the respective right or left derivative is infinite. For proofs of these 
statements in a general context we direct the reader to [|27l . 

In Proposition l3.4l we showed that a G (0, 00) belongs to t~^ (7) if and only if — log a 
is a subgradient of S at 7. Since r is nondecreasing and r(0) = 0, it follows from this 
result together with the preceding analysis that 

x-\0)= [o,e-^^(°)l, 
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t-Vil 



e~'^'i^^\+oo 



and for each 7 G (0, 1), 



e-^'A-y) (,-s'Ai) 



where e~°° is understood as zero and 6"*"°° is understood as +00. Our task in this section, 
then, is to compute these left and right derivatives explicitly in the case of rational 7 
(showing in the process that they are finite at and 1 if and only if the appropriate matrix 
is diagonalisable) and then show that for 76 (0, 1) fl Q the left and right derivatives of S 
at 7 cannot be equal to one another. 

To begin the proof we treat those statements concerned with t~^(0) and t^^(l). Let 
us suppose first that Aq is diagonalisable. Since Aq is non-negative it has an eigenvalue 
equal to its spectral radius (see e.g. [15. Theorem 8.3.1]), and since by Definition 12.21 its 
trace is positive, the remaining eigenvalue lies in the interval {— p{Aq) , p{Aq)\. It follows 
easily that the limit Pq := lim„_>.oo p(Ao)""Ag exists. Using Proposition 13 .4r iii) together 
with the fact that the word 0"1 is cyclically balanced, we may now calculate the right 
derivative of 5 at as: 

'^ ' n^oo l/(n+l) 

= lim (n + 1) (^ logp(A^Ai) - logp(Ao)) 
= lim logp(A[}Ai) - (n + l)logp(Ao) 

It follows that t^^(O) = [0,p(y4o)/p(Po^i)] as claimed in the statement of Theorem l2.3r ii). 
Let us now suppose instead that Aq is not diagonalisable. In this case Aq has a repeated 
eigenvalue equal to its spectral radius and has nontrivial Jordan form. It follows that 
lim„_^oo ||p(^o)~"^oll — +00. Let 5 > be the smallest entry of the matrix AqAi, 
which is positive by Definition I2.2r ii). and for each n > 2 let m„ be the largest entry 
of the non-negative matrix p(Ao)~"~^Aq^^ Clearly we have lim„_j.oo ^n = +00. Since 
Aq"^ and AqAi are both non-negative matrices it follows easily that 

2p (p(Ao)-"-M^Ai) > tr (p(Ao)-"-M^Ai) > 6mn 

for each n > 2, and hence 

S;(0) = lim 5(l/(" + D) - 3(0) ^ j.^ (^_A'^A,) = +0=. 

"^^ ^ n^oo l/{n+l) n^oo "^ ^ \p{Ao)''+^ " J 

In this case we therefore have c^^(O) = {0} as claimed. The proof of the statement 
concerning t~^(l) and the matrix Ai is almost identical, and we omit it for the sake of 
brevity. 
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Let us now move on to the case of 7 G (0, 1) fl Q. Fix p/q E (0, 1) for the remainder 
of the proof, and let {u, v) be a standard pair such that ^(ut>) = p/q, which exists by 
Proposition O Define Bi := A{u), B2 := A{v), and A := B1B2. By Definition Qii) 
the matrix A is positive, and hence by the Perron-Frobenius theorem A has two distinct 
eigenvalues. It follows that the limit 

P := lim p(A)-"A" 



n—^oo 



exists, and is the matrix corresponding to the unique projection whose image is the leading 
eigenspace of A and whose kernel is the non-leading eigenspace of A. In particular, P is 
of rank one. 

By Proposition 14.21 the word uvuv is cyclically balanced and the word u^v'^ is not: 
since both words have slope p/q and length 2q, it follows from Proposition 13 .4r iii) that 

p{B^B2y = piA{uvuv)) = e^'^^^P/") > piAiuh^)) = p{BlBl). (5.1) 

Using Proposition 13.40 1) we in particular have q{q) = e^^^l'^^a^l'^ = p{Aa{uv)y^'' for 
every a E ^^^(p/q) as claimed in Theorem 12.30 1). Let pi := \u\i, p2 := |f |i, qi := \u\, 
q2 := |f |. Using Proposition 13.40 11) together with the fact that the words (uv)"'v and 
u{uvY are cyclically balanced, we have for each n > 1, 

S ('^^P±^\ = ^^\og p{A{u{uvr)) = -J^ log p{B^ A"), 
\nq + qij nq + qi nq + qi 

S (!^l±^) =^^logp(A{iuvrv)) = ^^\ogpiA-B2). 
\nq + q2 J nq + q2 nq + q2 

Since 1^1 ■ It'll — |m|i ■ l^;! = 1 we have ^(-u) < ^(f ), and therefore q((uv)'^u) < (,{uv) < 
<^((f u)"f ) for all n > 0. We may therefore compute the left derivative of S at p/q as 

,p\ S{^-S{^^ liogp{A) -^-\og p{B,A-) 



p _ np+pi ^^^ p _ np+pi 

q nq+qi q nq+qi 

{nq + qi) log p{A) - q\ogp{BiA'^) 
= lim 

n-s>oo pq^ — PlQ 

where we have used the identity 

PI1 ^ PiQ = {\u\i + \v\i)\u\ — Iwlidwl + |t>|) = |-u| ■ |t>|i — ImIi ■ |t>| = 1. 
A similar calculation for the right derivative yields 



^^'''='°H7w 
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Combining this with the observations at the start of this section we obtain the explicit 
formula 



,-i,P 



.9 
asserted by Theorem l2.3[ 



^;(?) p-^.'(?)' 



[ p{A)i^ ' piPB,)'! 



(5.2) 



5+v^ 


V5 


10 


5 


V5 


5-V5 


5 


10 



Example 5.1. Let Aq, Ax be as in (11.11) and let p/g = 1/2; here li = 0, f = 1, whence 

Pi = 0, gi = l,p2 = 1, g2 = 1- Therefore, 

whence by (15.21) . r^^(l/2) = [4/5, 5/4]. This result was previously derived in [|30l by a 
different, geometric method. (See also Table 12. H and Figure [IJ) 

To finish the proof of Theorem l2.3r ii) we need to show that the interior of t~^(p/g) is 
nonempty for any p, q. Suppose it is empty; then 

p{BrPYp{PB2Y = p^AY'^"^ = p{AY. (5.3) 

For the remainder of the proof we shall assume that this relation holds, and thereby derive 
a contradiction. 

So, suppose p{BiP)p{PB2) = piA). It follows from the definition of P that p(A) = 
p(PAP). Since the spectral radius of a product of matrices is invariant under cyclic per- 
mutations of that product, we have p{PBi) = p(BiP) and p(Pi?2) = p{B2P), and it 
follows that p(PBi)p(B2P) = p{A) = p{PAP). Since A is invertible it has nonzero 
spectral radius, and therefore p{PBi) and p(i?2-P) are positive. It follows from the def- 
inition of P that P is a non-negative matrix, and hence PPi, B2P and PAP are all 
non-negative. In particular, each of these three matrices has an eigenvalue equal to its 
spectral radius. On the other hand since P has rank one, each of the matrices PPi, B2P 
and PAP has determinant zero, and hence has one eigenvalue equal to zero. We conclude 
that (trPPi)(trP2^) = tr PAP. 

Since the two-dimensional matrix PPi has one positive eigenvalue and one eigenvalue 
equal to zero, it is diagonalizable. Define Ai := p(PPi) and A2 := p(P2-P), and choose 
an invertible matrix U such that 



where a, 6, c and d are real numbers. We have A2 = tr B2P = tr UB2PU^^ = a + d 
and A1A2 = (trPPi)(trP2P) = tr PAP = ti PB1B2P = ti UPBiU-^UB2PU~^ = 
Xia. It follows that d = and a = A2. On the other hand, since det P = we have 



20 Ian D. Morris, Nikita Sidorov 

det UB2PU'^ = and therefore ad — he = Q. We deduce that he = 0, and therefore at 
least one of h and c is zero. 

We claim that there exists a nonzero vector w G M^ which is an eigenvector of both Bi 
and i?2. We consider separately the case 6 = and the case c = 0. 

If 6 = 0, put u = f/^HO, 1)^, and note that 



whence 

PBioo = U- 

and 

B2PUJ = u- 







Ai 0\ /0\ ^ /O 
Oj{lJ {0 

': t) (?) - C 



i.e., PBicu = B2PUJ = 0. We claim that to is an eigenvector of Bi and of B2. Since B2 
is invertible we deduce from the relation B2PU1 = that to belongs to the kernel of P. 
Since Bi is invertible and PBi = 0, we deduce that Biu also lies in the kernel of P, 
and therefore u is an eigenvector of Bi. Finally, since by definition the kernel of P is one 
of the eigenspaces of A = B1B2, the vector cu is also an eigenvector of B1B2, and since 
B^^cu is proportional to to we conclude that B2UJ is proportional to tu as required. 

The case c = is similar. Define to := U^^{1, 0)^ so that PBitu = Xico and B2PU1 = 
X2U!. Since to = Xi^PBicu the vector cu belongs to the image of P, and is therefore fixed 
by P since P is a projection. It follows from this and the relation B2PUJ = A2W that 
B2OJ = A2W, so that cu is an eigenvector of P2- On the other hand, the image of P is an 
eigenspace of A and therefore cu is an eigenvector of A = B1B2. Since B2CU = \2cu we 
deduce from this that u is also an eigenvector of Pi as required. This proves the claim. 

We may now derive the desired contradiction. Let w be a common eigenvector of the 
matrices Pi and P2, and let us write PjW = ^iU for i = 1,2. We have BlB^cu = ^i^fw = 
{B1B2YU} so that ^1^1 is an eigenvalue of both Pi P| and {B1B2Y ■ It follows immediately 
that p{BlBl) = max{|e?ell, I (det PiP2)2^r^^2"^l} = p{BiB2f, contradicting (EH)- We 
conclude that the relation (15.31 ) cannot hold, and the proof of Theorem |2.3r ii) is complete. 

Corollary 5.2. The function S : [0,1] -^ Wis strictly concave. 

Proof. We know that S is concave; if it were not strictly concave, there would be an 
intervalJ = (71,72) suchthat^jj would be affine, i.e., S"|j would exist(and be constant). 
This contradicts the fact that S'i^j) does not exist for any rational 7. D 

Remark 5.3. One can show that for the matrices Pi and P2 we have 

p(PlP2) > p(Pl)p(P2), 

which is essentially equivalent to Corollary 15.21 We leave this as an exercise for the inter- 
ested reader. 
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6 Preimages of irrational points: preliminaries 

In this section we apply a dynamical estimate to prove an inequality dealing with the 
subgradients of the function S. We begin with the following lemma, which allows us to 
choose a family of norms on M? which is adapted to the study of the family Aq,. The proof 
is identical to that of [fT3l Lemma 3.3]. Note that the proof requires that Aq and Ai do not 
have a common invariant subspace, as was stipulated in Definition 12. 2r i). 

Lemma 6.1. There exists a family of norms { || ■ ||q, : a G (0, oo)} on M^ with the following 

properties. For every v G M^, a G (0, oo) and i = 0,1 we have WA^'vUa < f?(a)||f ||q,. 
If K G (0, oo) is compact, then there is a constant C > 1 depending on K such that 
\\v\\ < C\\v\\a < C^\\v\\forallv G M^ and a G K. 

For the remainder of this section we fix a family of norms || ■ ||q, with the above prop- 
erties. We also make use of the following elementary result: 

Lemma 6.2. Let B be a 2 x 2 invertible real matrix, and let \\ ■ \\ be the Euclidean norm. 
Then there exists a rank one matrix Q such that \\B — Q\\ = \ det i?|/||i?||. 

Proof. We apply the existence of a singular value decomposition for B. Let us choose 
unitary matrices U, V and a non-negative diagonal matrix D such that B = UDV^. Since 
U and V are isometrics with respect to the Euclidean norm we have \\D\\ = \\B\\, and 
since | det f/| = | det V\ = 1 we have | det D\ = \ det i?|. It follows that the nonzero 
entries of D are ||i?|| and | det i?|/||S||. If P is a matrix which has the entry ||i?|| in the 
same position as for D, with all of its other entries being zero, then clearly P has rank 
one and \\D - P\\ = | det 5|/||5||. Now let Q := UPV^. D 

We now require a dynamical result which describes the dependence of the eigenvec- 
tors of certain products Aa{u) on the structure of the words u. The following result is 
similar in spirit to [|23l Theorem 2.2], but has the additional property that the modulus of 
continuity of the vector-valued function depends on a in a controllable manner. The re- 
striction to rational 7 serves only to simplify the proof: by working instead with two-sided 
Sturmian sequences indexed over Z, this condition could be removed. The rational case 
being sufficient for our argument, we ignore the more general statement. 

Proposition 6.3. Let L C (0, 1) be compact, and let 7 = t{a) G L fl Q. Then there exist 
constants 6 G (0, 1) and K > 1 depending only on L and a function : X^ — )■ M? 
such that the following properties hold. For each x G X^, t)(x) is non-negative and 
satisfies \\v{x)\\a = 1 and Aa{x,n)\3{x) = g{a)"'\3(T"'x) for all n > 1. Ifx,y G X^ 
with d{x,y) < 2'^ then \\\3{T^x) - t)(T'=y)|U < KO''. 
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Proof. Let us write 7 = p/q in least terms. Consider any x E Xj. since < p/q < 1 the 
matrix Aaix, q) is a mixed product of Aq and Ai, and so is positive by Definition 12 ■2r ii). 
By the Perron-Frobenius theorem it follows that Aa{x., q) has a unique positive eigenvec- 
tor with associated eigenvalue equal to p{Aa{x, q)) = limfc_i.oo ||-4.q,(x, kq)\\^/^ = g{a)''. 
Let i){x) be a positive vector belonging to this eigenspace such that \\V){x) \\a = 1. 

Let us fix X G Xj and show that Aa{x, n))o{x) = g(a)"'t>{T^x) for all ra > lAfn>q 
then 

Aa{x, n)X){x) = AaiT'^x, n — q)Aa{x, q)'o{x) 
= g{ayAa{T'^x,n-q)X){x) 
= g{ayAa{x,n — g)t)(x), 

and by iterating this identity we may reduce to the case where 1 < n < g. In this case we 
have 

g{ay = \\Aa{x,q)v{x)\\a = \\AaiT^x,q — n)Aaix, n)v{x)\\a 

< \\Aa{T''x,q - n)\mAa{x,n)v{x)\\a 
<Qiay-''\\A^{x,nMx)\\^<g{ay 

and it follows that ||^Q(a;, n)t>{x) \\a = g{a)^. Moreover we have 

AaiT'^x, q)Aa{x, n)\3{x) = Aa{T'^x, n)Aa{x, g)o(x) 

= g{ayAaiT'^x,n)i){x) = g{ayAaix,n)v{x). 

Thus Aa{x, n)X){x) is a positive eigenvector of ^a(T"a;, q) with corresponding eigenvalue 
g{ay and with norm equal to g{a)"'. Using the uniqueness of the leading eigenspace in the 
Perron-Frobenius theorem we conclude that g{a)~'^Aa{x, n)t)(x) = t)(T"x) as claimed. 
Since r is continuous and monotone with r(0) = and limQ,_j.oo ^{ce) = 1, the set 
r^^(L) is a closed subset of (0, 00) which is bounded away from and cxd, hence compact. 
For every a E x:^^{L) we have m.ax{p{Ao), p{aAi)} < g{a) by Definition I2.2f v). since 
x{a) ^ {0, 1}. Let us define 

(\detAo\ |det(aAi)n max{p(Ao)^pMl)2} 
e := sup max <^ — --— , -— — )■ < sup —-^ < 1, 

and let x,y E X^ with d{x, y) < 2"^ We will show that ||d(T'=x) - X){T^y)\\a < QC^O^, 
where C > 1 is the constant provided by Lemma 16.11 with respect to the compact set 
x~^{L)Ai?,C^e^ > 1 then clearly \\^{x) -o(y)|U < 2 < GC^^'^, so we shall assume for 
the remainder of the proof that SC^O'' < 1. Since 

\det{g{a)-'^Aa{x,k))\ ^ C'O^ _ ^,^, 



\g{a) ^Aa{,x,k)\\ \\g{,a) ^Aa{,x,k)\\o, 
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it follows from Lemma \62\ that there exists a rank one matrix Q E M2(M) such that 

\\Q{a)-''Aa{x,k)-Q\\ < C^ 9''. Since d{x,y) < 2-^' we have ^„(x, /c) =Aa{y,k) and 
therefore also Wgiay^Aaiy, k) - Q\\ < C^O^. Clearly, 

\\Q{a)-^A^{x, kMx) - Qo(a;)|U < C^^' < ^ < ^' 

and therefore in particular |1 - ||(5o(2;)||a| < 1/3 so that 2/3 < ||Qo(x)||q, < 4/3. 
By identical reasoning we have 2/3 < ||(5t3(2/)||a < 4/3. Now, since the image of Q 
is one-dimensional, there exists A G M such that Qti{x) = \Q'o{y). We have |A| = 
\\Qxi{x)\\a/\\Q^{y)\\a < 2 and therefore 

||0(T'=x) - Ad(T^)|U = \\Q{a)~^A^{x, fc)t)(x) - Af?(a)-^A(i/, fc)o(y)|U 

< \\Q{ay^Aa{x,k)x){x) - Q\3{x)\\o, 
+ \\\Q\3{y) - \Q{a)-^Ac,{y, k)x>{y)\\a 

< (1 + |A|)C'^'^<3C'^^ 

We now claim that the real number A is positive. Suppose that it is negative: since X){T^x) 
dcad.X){T^y) are non-negative, we must have \\\3{T^x) — \'o{T^y)\\ > ||t)(T^x)|| and there- 
fore 

||o(r'=x) - Ao(T'=x)|U > C-^\\v{T''x) - AD(T'=y)|| 

> C-1d(T^x)|| > C-2||o(T^x)|U = C-2, 

which contradicts our assumption that 3(7^6*^ < 1. We conclude that A must be positive, 
and therefore 

|1-A| = |||D(T^a;)|| - ||AD(T'=y)|| I < ||d(T'=x) - AD(T'=y)|| < SC'^O''. 

It follows easily that p{T^x) - D(T'=i/)||a < QC^O'' < QC^9^. The proof of the proposi- 
tion is complete. D 

Finally, we make use of the following simple result from matrix analysis, which we 
adapt from [9, Lemma 2]. 

Lemma 6.4. Let ||| ■ ||| be a norm on M^, let B be a 2x2 matrix with \B\ < 1, and suppose 
that V eM? with \v\ = 1. Let C > 1 be any constant such that C~"^||m|| < |||m||| < C||m|| 



for all u e W. Then 1 - 2C^^l\Bv - v\l < p{B) < 1. 

Proof. If Ml and M2 are a pair of 2 x 2 real matrices, /i is an eigenvalue of M2, and Ai, A2 
are the eigenvalues of Mi, then the bound 

min{|Ai - /i|, IA2 - /i|} < v/(l|Mi|| + ||M2||)(||Mi-M2||) 
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is well-known, see for example [2, §VIII]. Define Mi := B and M2 := B + ||f ||~^(f — 
Bv)v^; we may then estimate 

||Mi - M2II = ||t;|r^.||(t^ - Bv)v^\\ < \\v\\-\\\Bv - v\\ < C^jBv - vj 

and 

||Mi|| + IIM2II < 2||Mi|| + IIM1-M2II < 2C^l\B\l + C^llBv - v\l <AC^. 

Since M2V = t>, 1 is an eigenvalue of M2, and it follows that B has an eigenvalue A such 



that |A — 1| < 2(7^ A/||5t7^^t7| . The result follows in view of the elementary inequality 
|A| < p(fi) < |||fi||| = 1. D 

The following key estimate forms the core of the proof of Theorem I2.3r iii): 

Lemma 6.5. Let L C (0, 1) be a compact interval, let 7 = Pn/Qn = [^i, • • • , c^n] ^ L 
where n,an > 1, and choose any a G (0, 00) such that r(a) = 7. Let Pn~i/<ln~i = 
[ai,. . .,a„_i]. Then 

Nn+l fS{l) - S{pn-l/qn-l) 



V 1-Pn-l/qn-l J 

where K > 1 and 6 & (0, 1) are constants depending only on L. 

Proof. From the classical theory of continued fractions we have 7 < Pn-i/qn-i if n is 
even, and the reverse inequality holds if n is odd. By Proposition 13 .41 we know that S is 
concave and (— log a) is a subgradient of 5 at 7. Since the average gradient of a concave 
function on a closed interval is bounded below by every subgradient at the right endpoint, 
and bounded above by any subgradient at the left endpoint, we immediately deduce the 
inequality 

^n+l fS{-f)-S{pn-l/qn-l) 



V 7 - Pn-l/qn-l 

In proving the remainder of the lemma we will assume that q-n > rn for some constant 
m> 1 to be determined below. Indeed, given any such m it is clear that L contains only 
finitely many rational numbers Pn/qn with denominator less than m, and so by adjusting 
the constant if > 1 if necessary, the full strength of the lemma follows from this special 
case. Since the hypotheses of Proposition 14.11 are satisfied, we may fix an integer k > 
|g„ and a point x e X^ such that d{x, T'^"-^x) < 2"^ and q^\ log p{A{T^x, qn-i)) = 
S{pn-i/qn~i)- Let y := T'^x. By Proposition l6.3l we have 

for some constants K > 1 and 9 G (0, 1) depending on L. Combining this with Lemma [64l 
we obtain 

i-2C^VKe-^'^" < ^(a)-^"-v(A(y,gn-i)) < 1. 
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where C > 1 is the constant assigned by Lemma [611 to the compact set t^^(L) C (0, oo). 
Since x{a) = 7 we have g{a) = e^^"''>a'^, and therefore 

Let m be an integer which is large enough that 1 — 2C'^\/K6'^^^ > e~^. By considering 
only those cases where g„ > m, we may by taking logarithms obtain 

-2C^y/Ke-e''" < q„,-iS{pn-i/qn-i) - qn-iS{l) + [Pn-i - Qn-il) log a < 0, 
and by a slight rearrangement, 

< 5(7) - S{pr.-ilqn-i) +{1- ^-^\ log a < '^-^^^^eh-. 

V (In-lJ qn-1 



,n+l/., Pn~A i-^r^\^ . „^ ^ 1 



Now, since 

V Qn-1 J Qn-lQn Qn-lQn 

we may derive the inequality 

< (-1)"+^ p(7)-^(Pn-i/9n-i) ^^^g \ ^ 2qr,C'^/Keh-, 

V 7-Pn-l/gn-l / 

which completes the proof. D 

The following interesting result may be derived from Lemma [631 

Corollary 6.6. Let L C (0, 1) be compact. Then there exist constants C > 1, ^ G (0, 1) 
depending on L such that for all p/q G L with p and q coprime, the length of the interval 
t~^{p/q) is bounded by CqO'^. 

Proof. By enlarging the constant C if necessary it is sufficient to consider rationals p/q 
which are not of the form 1/k for an integer k >\, since L can contain only finitely many 
rationals of this form. Let p/q E L with p and q coprime and p > 1. These assumptions 
allow us to find an integer n > 2 and integers ai, . . . , a„ with a„ > 1 such that p/q = 
[ai, . . . , a„]. If «! and 02 are the endpoints of the closed interval t~^{p/q), then we may 
apply Lemma [631 twice with a = ai, 0:2 to see that | log a2 — log ai | < 2Kq9'^. The result 
follows. D 

Recall that if 7 = [oi, 02, . . .] G (0, 1) \ Q and the sequence (a„) is bounded, then 
there exists a constant 5 > such that I7 — p/q\ > 5q^'^ for all g G N and p E Z (see 
for example ([TSl ). In particular, for all sufficiently large k the relation I7 — p/q\ < \/(^ 
is impossible for integers p E TL and q > 2. The proof of the following lemma is thus 
identical to the proof of [13, Lemma 8.3]: 

Lemma 6.7. Let 7 = [ai, 02, . . .] G (0, 1) \ Q, and suppose that a„ = 1/or all sufficiently 
large n. Then t^^(7) is a singleton set. 
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7 Preimages of irrational points: proof of Theorem. 

Let us define Z := r~^((0, 1) \ Q), and partition Z into two subsets as follows. We define 
Zi to be the set of alia E Z such that the infinite continued fraction expansion [ai,a2, ■ ■ ■] 
of the irrational number r(a) G (0, 1) satisfies a^ = 1 for all but finitely many k. If 
7 = [ai, 02, . . .] is an irrational number of this type, then by Lemma [677] the set t^^(7) is 
a singleton set. It follows that Zi is countable, and hence has zero Hausdorff dimension. 
Let us now define Zq := Z \ Zi. Since the Hausdorff dimension of a countable union 
of sets is equal to the supremum of their individual Hausdorff dimensions, to prove that 
diiRHiZ) = as claimed it is sufficient (and indeed necessary) to show that dimj|/(Zo) = 
0. Moreover, it is sufficient to show that for some sequence of sets Lk whose union covers 
(0, 1), each of the sets Zq fl t^^(Lfc) has Hausdorff dimension equal to zero. For the 
remainder of the proof, we fix a set L of the form [i, 1 — i] with the aim of showing that 
the set tr^^(L) fl Zq has Hausdorff dimension zero. 

Given natural numbers n, Oi, 02, . . . , a„, let T(^ai,...,a„) denote the half-open interval 
with endpoints [oi, . . . , an] and [oi, . . . , 1 + a„] which excludes the former endpoint but 
includes the latter. This interval consists precisely of those elements of (0, 1) which admit 
a continued fraction expansion whose first n entries are ai, . . . , a„ respectively, and whose 
length is at least n + I. Given natural numbers n > I and ai, . . . , a„, let us define 

I(ai,...,a„) ■= ^' {L n T(^ai,...,a„)) ■ 

Note that by our choice of L, if n > 2 then the set I(ai,...,a„) is either empty or is equal to 
all of x^^{T(^ai,...,an))- For s^'^h A^ > 2 let us define Un to be the set of all 'I(ai,...,an) such 
that a„ > 1, ra > A^, and a^ = 1 for all k such that N < k < n. The reader may easily 
verify that T(^au-,an) ^ ^{bi,...,bm} = when the vectors (ai, . . . , a„) and (61, ... , fe^) are 
distinct, and furthermore, 

t-\L)nZoC U X^a,_a^) (7.1) 

(ai,...,a,i)gl/]v 

for every N > 2. We make the following key claim: if n > 2 is an integer, then for each 
n-tuple of natural numbers (ai, . . . , a„) G N" such that a„ > 1 there holds the inequality 

diamX(„,,...,,„) < KqJ'^", (7.2) 

where Pn/qn '■= [«i, ■ ■ ■ ,o.n\ in least terms, and K > 1 and 9 G (0, 1) are constants 
depending only on L. 

Let us prove this claim. Fix an integer n > 2 and suppose that T(^ai,...,an) is nonempty 
with a„ > L Let ai and 0:2 be respectively the infimum and the supremum ofI(^ai,...,an), 
and let 7j = x{ai) for i = 1, 2. If n is odd then we have 71 = [ai, . . . , 1 + a„] and 
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72 = [oi, • • • 5 o„], and if n is even then 71 = [oi, . . . , a„] and 72 = [ai, . . . , 1 + a„]. 
Define also pn-i/qn-i '■= [^i, . . . , a„_i] = limj_>.oo[ai, • • • , On-i, "^l- Our objective is to 
bound the difference ^2 — «!• 

We consider first the case in which n is odd, in which case Pn-i/qn~i < 7i- Recall that 
for a concave function defined on an interval [a,b], the average gradient in the interval is 
greater than the value of any subgradient at b, and less than the value of any subgradient 
at a. Since (— log 02) is a subgradient of S at 72, and (— logai) is a subgradient of S at 
7i, it follows that 

'5(72) -'S'(7i) > (72-71) (-log "2) 

and 

S'(7i) - S'(p„_i/g„_i) > (71 -p„_i/g„_i)(-logai). 

Adding these two inequalities together, we obtain 

(72 -71) (-log "2) + (71 -Pn-l/gn-l)(-logai) < S'(72) - S{pn-l/qn~l) 

and therefore 

(71 -Pn-i/gn-i)(loga2 - log«i) <5'(72) - 5'(p„_i/g„_i) 

+ (72-Pn-l/gn-l)loga2. 

Since L is a compact subinterval of (0, 1), t^^(L) is a compact subinterval of (0, 00), so 
there is a constant C > depending on L such that | log x — logy\ > C^^ \x — y\ for every 
x,y e L. Hence 

C^^{a2 — ai) < loga2 — logai 

^ l' l2-Pn-l/qn^l \ / g(72) - S{pn-l/ Qn-l) , Jq ^ A 
~ \ll-Pn~l/qn-l) \ l2-Pn-l/qn-l V' 

Let Pn-2/qn-2 = [^1, • • • , a„_2] iu least terms. Since 

_ (1 + an)Pn-l +Pn^2 _ dnPn-l + Pn-2 

7l fT \ \ I ' '•^ I ' 

U + (In)qn-l + qn~2 CtnQ'n-l + Q'n-2 

and qn^iPn-2 - qn~2Pn-i = 1 it follows that 

72-Pn-l/gn-l _ l/(QngLl + gn-ign-2) 

7i - Pn-i/qn-1 1/((1 + an)ql-i + g„-ig„-2) 



1 + ar, 



qn-2 

q-n-i 



< ^ < 3. 

On 
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Applying Lemma [631 we obtain 

^2 - «! < 36 (-1)"^' h log ^2 

V l2-Pn~l/qn-l J 

< SCKqJ'^" 

as required, which completes the proof of the claim in the case where n is odd. 

We now consider the case in which n is even. In this case we have Pn-i/Qn-i > 72- By 
comparing subgradients in a similar manner to the odd case we arrive at the inequalities 

^(72) -^(71) < (72-7i)(-logai), 

S{Pn-l/qn-l) - 5'(72) < (Pn-l/qn-l " 72)(-loga2)- 

Adding these two inequalities yields 

SiPn-l/qn-l) - S{-fi) < (Pn-l/qn-l - I2) {- '^Og a2) + (72 -7l)(-logai) 

and therefore 

S{Pn-l/qn-l) -5'(7i) + {Pn^l/qn-l -7l)logttl < (72 " Pn-l/^n-l) (log ^2 -logai). 

Dividing by the negative real number 72 — Pn-i/qn-i we obtain 



1 1 ^ f-fi-pn~iqn~l\fS{-fi)-S{pn-iqn--l) ,, 

log aa - log ai < - h log ai 

\l2-Pn-l/qn-lJ \ ll-Pn-l/qn-1 

^3(_,).^ip(7i)-g(Pn-i/gn-i)^^^g^^ 

V ll-Pn-l/qn-1 

and it follows using Lemma [631 that 02 — ai < 3CKqn9'^" as before. This completes the 
proof of the claim. 

We may now show directly that t~^(L) fl Zq has Hausdorff dimension zero. We re- 
call the definition of the Hausdorff dimension of a set F C R. For each A > 0, the 
A-dimensional Hausdorff outer measure of the set Y is defined to be the quantity 

hm inf < NJ(diam U)^ '■ Y C \\ U and sup diam U < 6 > , 

where each W is a collection of subsets of R. The Hausdorff dimension of the set Y is then 
defined to be the infimum of the set of all A > such that the A-dimensional Hausdorff 
outer measure of Y is zero, or equivalently the infimum of the set of all A > for which 
this value is finite. 
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Let A G (0, 1], and choose any 5 > 0. We saw in (|7.1I) that the union of the elements of 
Un contains tr^^(L) fl Zq for every A^ > 2. It follows from (17.21) that if A^ is large enough 
then every element of Z^at has diameter less than 6. For any such A^ we have 

J2 (diamX(„,...,„„))' < Yl i^^lnO'")' 

(ai,...,a„)eWjv (ai,...,an)6Wjv 

< Yl (^^^')' 

p/ge(Qn{o,i) 

9=2 p=l 

< KYq^e^'i = ^-— — J — ^, 



9=1 



and since this bound is independent of 6, we conclude that the A-dimensional Hausdorff 
measure of r^^(L) flZo is finite. Since A may be chosen arbitrarily close to 0, we conclude 
that dimj|/(r"^(L) fl Zq) = as required. The proof of Theorem |23l iii) is complete. 

8 Explicit formulae 

In this section we prove Theorem 12.51 and present some bounds which can be used for 
practical computation of t~^(7) in the special case of the matrices defined by (11.11) when 
7 is not too well approximated by rationals. In [13] we proved the following result (the 
indexing of the sequences in the statement of Theorem 18.11 has been adjusted so as to 
agree with the conventions used elsewhere in this section): 

Theorem 8.1. Let (T„)5^g denote the sequence of integers defined by r_2 := 1, r_i, tq := 
2, and 

^n+l := TnTn-l - T,a-2forall Tl > 0, (8.1) 

and let {Fn)'^=Q denote the sequence of Fibonacci numbers, defined by Fq := 1, Fi := 1 
and Fn+i := F„ + -F„_i for all n > 1. For each a > let Aq, be the pair of matrices 
defined by (|l.ll) . and define a real number a^, G (0, 1] by 

Then this infinite product converges unconditionally, and Aq, does not have the finiteness 
property. The numerical value of the constant a^, is 

a* ~ 0.74932654633036755794396194809. . . 
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Here a* is in fact the unique positive real number such that x{a^:) = (3 — y/b) /2. This 
particular constant was studied because % := (3 — \/5)/2 has a particularly simple con- 
tinued fraction expansion: we have 7^, = [2, 1, 1,1,1,... ], which is the simplest possible 
expansion of an element of (0, 1/2) \ Q. 

Now that Theorem 12. 3| has been proved, the proof of Theorem [23] may be obtained in 
a manner essentially similar to the proof of Theorem 1 8. II 

Proof of Theorem |23] Let 7, (g„), (s„) and (p„) be as in the statement of the theorem, and 
let r be as in Definition l2.2[ An inductive argument as used in §4 shows that | s^ 1 1 = Pn and 
\sn\ = gnforeveryra > 1. In particular we have ^(s„) = p„/5'n for alln > 1 and therefore 
q~^ log Pn = S{pn/qn) for every positive integer n by Proposition l3.4[ By Theorem l2.3r i) 
the function r is continuous, monotone non-decreasing, and satisfies r((0, +00)) 3 (0, 1). 
In particular, 1:^^(7) is nonempty, and is either a point or a closed interval. A consequence 
of Theorem I2.3r iii) is that t~^(7) has empty interior, and we conclude that there is a 
unique point a^ G (0, 00) which satisfies t(a^) = 7. It follows via Proposition 13.41 that 
— log a^ G M is the unique subderivative of S at 7, and hence S is differentiable at 7 with 
5" (7) = — log a^. We may therefore calculate 



S"(7) = lim 



S {"^^^ - 5 f £ii 



<ln + l <ln 

^— log Pn+l log Pn 

n-XL P2i±l_Pn 

Qn + l Qn 

qnlogpn+l-Qn+l'^Ogpn 

= lim 

n^oo QnPn+l " Qn+lPn 

= lim (-l)"(g„logp„+i - g„+ilogp„), 

the existence of all of these limits being guaranteed by the differentiability of S* at 7. By 
rearranging we obtain 

Pn+1 

as claimed. To derive the product expression for a^ let us define 

(-1)" 



ttr 



/iQn + l 

rn 
Pn+1 



for each n> —\, and observe that 

"n / Hn Hn-1 \ Hn Hn-1 \ Pn Pn-1 



"n-1 \Pn+lPn J \ Pn+1 J \ Pn+1 
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for each n >0, where we have used the relation g„+i = a„+ig„ + Qn-i- We also have 

pV PiAof 1 



a_i 



p'i\ piA.y piA, 



Hence 



a^ = lim a]\r = lim a_i TT 






n=0 



lim , _, , TT 



-■- TT / Pn Hn-l 



N^oo p{Ai) A^A \^ p^^^ 



•^ / ^a„+i^ \ (-l)"gr, 



-■- TT I "n Fn-1 



as claimed. The proof is complete. D 

For the remainder of the section we let { A^ : a > 0} be the specific family of matrices 
defined by (11.11) . In this case we have p{Ai) = 1, which means that the term 1/ p(Ai) 
may be removed from the infinite product formula in Theorem 12.51 Let 7 G (0, 1) with 
infinite continued fraction expansion given by 7 = [ai, 02, a^, . . . ], and let Pn/qn be the 
nth convergent of 7. In view of the identity Aq = Ai, by replacing 7 with 1 — 7 and 
a with 1/a if necessary, we will assume without loss of generality that a E (0, 1) and 
7 G (0, 1/2), which is equivalent to ai > 2. 

Let us consider the sequence of words specified by 7 given by s_i = l,so = 0, 
si := Sq^^-'^S-i and s„+i = s^"+is„_i for all n > L Since s„ prefixes s„+i for every n, it 
follows that Sn prefixes Sk for every k > n. Since furthermore the lengths |s„| = g„ tend 
to infinity, it follows that there is a unique infinite word s^o G ^2 which is prefixed by 
every s„. In particular this word is balanced, and it is recurrent: for each n > the prefix 
s„_i occurs in at least two distinct locations in the prefix s„+i, hence at least four distinct 
locations in the prefix s„_|_3, and so forth, so that every sub word of Sqo recurs in infinitely 
many positions. 

Since <i{sn) = Pn/ln, we have ';(s„) -> 7 as n — )■ 00 and using Theorem 12.11 it 
follows that Soo G X^. For 7 = 7* the word Soo is none other than the Fibonacci word 
010010101001 . . . , which is the fixed point of the substitution 0^01,1^0 (see IfTOj 
|22i). 

Define i?„ = A{sn) for each n> -1. We have i?_i = Ai,Bo = Aq, Bi = A^^^'^Ai, 
and 

S„+i = S^"+^S„_i, n > 1. (8.3) 

Put Tn = tr Bn and p„ = p{Bn) as before. Note that since det i?„ = 1 we have r„ = 
Pn + Pn^ and conversely pn = \{Tn + a/t^ — 4) . In particular r„ ~ p„ as n — )■ 00. 
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Subject to the above hypotheses we will prove the following rigorous estimate for the 
error in approximating a^ by a partial product: 

Proposition 8.2. Suppose there exists a constant L > and an integer no > 3 such that 

In < Lpn-ifor alln > uq. (8.4) 

Then for every N > uq there holds the inequality 



^ / .a„+i ^ \ (-l)"'2r, 



'---'Hn(^^) 



< ^. (8.5) 

Pn 



where Cq := 16(ai + l)(ai + 2) + 1. 

Remark 8.3. The assumption (18.41) is very weak. In particular, it holds for any non- 
Liouville 7 - see Lemma [877] below. 

In the special case where the continued fraction coefficients of 7 are bounded, Propo- 
sition [8]2] lends itself to particularly easy verification. We have: 

Corollary 8.4. Suppose there exist integers K > 2 and tiq > 3 and a constant L > 
such that the inequalities qno+i < -^Pno '^'^d supjofc '■ k >2 + uq} < K — 1 are satisfied, 
and such that the matrix B„^_i — K ■ I is non-negative, where I denotes the identity. Then 
(18.51) holds for every N > uq. 

Since the spectral radii p„ grow super-exponentially as a function of n (see Lemma [83] 
below), this allows very exact estimates to be made using relatively few terms. In order 
to prove the proposition and its corollary we require two lemmas. The following result is 
the technical core of the proof: 



Lemma 8.5. The inequality 

Pn+l 



Pn Pn-l 



< -^ (8.6) 

Pn~l 



holds for alln > 1, where Cq is as in Proposition \8.2\ In particular, p„+i ~ p^"+^p„_i as 
n — 7- 00. 

Proof. We first construct an auxiliary continued fraction as follows: 

/3 = [rfi, 1,^2, 1,4, 1,---], 

where di = ai — 1 and dk is the number of zeros between the A;th and [k + l)st unities in 
Soo for k > 2. For instance, for ai = 2 and a^ = 1 for A; > 2 (i.e., the Fibonacci word Sqo) 
we have f3 = [1, 1, 2, 1, 1, 1, 2, ... ]. We denote 

/5=[&1,&2,...]- 



A Devil 's staircase from joint spectral radii 33 



Note that since the number of consecutive zeroes in s^o is bounded by ai (see, e.g., 112211 '). 
we have hk <ai for all k. 

Let Mo = 1, and for each n > 1 let m„ denote the length of the word constructed from 
s„ by replacing every string of consecutive zeros with a single zero. That is, si = 0"^^^1, 
whence ui = 2; S2 = {0"'^~^l)°''^0, whence U2 = 2a2 + 1, etc. Define also 

p 

-^ = [bi,b2,...,bk]. 

Recall the following well known relation between matrix products involving powers of 
Aq, Ai, and continued fractions: 

am Ad-m-l /lai _ I -L U \ / -'- '^m-l\ / -L U^ _ ( Pm Pm-1 



Aam A^m-l _ _ _ Aa 

if m is odd, and 

Aam Aam-1 A a 

^1 ^0 ■ ■ ■ /li 



am AO-m-l /lai _ / J- 0'm\ ( -L U^ /^i U^_ j Pm-1 Pn 



1 y yOm-i 1/ yai ly \Q'm-i Q'n 
if m is even (see, e.g., lEl). Hence we have 



Bn 



Un is odd 



(8.7) 



P P 

P P 

^Un — 1 ^Un 

Let us compute the eigenvectors for B.^ of the first type: 

t «::;) (t) - ^" (t) . <«■«) 

where A„ = p„ or p~^. Solving this system, we get a quadratic equation: 

^Un^n ~r W«n-1 J^^Unj^n ^M,i-1 ^ U. 

Dividing it by Qu^ , we obtain 
whence 

f en + %^ V^n - /3) = en f ^ - /?) + TT^ f ^ - /^V (^.10) 

Let from here on ^„ stand for the positive root of (18.91 ). From ( 18.81) it follows that Qu„^n + 

Qu,,-1 = Pn< Tn = Pu„ + Qu^-l, whcUCC ^„ < ^ < L 
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Since the bk are bounded, we have 

Qk . 1 



Qk-i 



bk + 



<bk + l <ai + l. 



(8.11) 



Jk+l 



Hence from (18.101 ) 



Kn-/3|<(ai + l) 



v«„ 



/3 



+ 



«n-l 



Wm„-1 



/3 



By (18.111 ). we have 



/3 



Hence 



g. 



< 



V«n*^«n + l 



< 



02 ' 



/3 



Pun-l 



V«n 



< 



^Un^ 



< 



Un — 1 



ai + l 



(8.12) 



i/3-e«i< 



(ai + l)(ai + 2) 



Since p„ < r„ = P„,^ + (5«„-i, we have p„ < 2(5„„, whence 

|/3-en|<%, 



(8.13) 



where 



Ci = 4(ai + l)(ai + 2). 

(In the case of even m„, we have r„ = -Pu„-i + Qu^ < '^Qu„, so (18.131) holds as well.) 
Let ^^ < denote the other solution of (fOt . Put 



We have 



Dr 



D-^BnDn 



1 1 



Pn 



Pn'J' 

We want to apply the change of coordinates given by D„ to the equation (18.31) and then 
obtain a relation for the traces. Since tr [D^^Bn+iDn] = tr -Bn+i = Tn+i = Pn+i + P„^+i> 
we will be only concerned with estimating tr (D^^Bn-iDn). 
Assume that Un-i is even; then 



B„-i 



P, 



Mn-l-l -* «n-l 



P„ 



(The case of odd n„_i is completely analogous.) We have 



Dn Bn-lDn 



1 



Sn 



Mn-l-1 -* Wn-1 



P„ 



Pn-1 - '"n-l • • 



Sn Sn 
1 1 



'n-1 
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where 

^n-l = ^ni^nQu„_i-l " -P«„_i-l) + (.nQu^^i - -P«„_i- (8.14) 

By (I8J21) and (I8l3]) . 



Pn Wwn-l 



+1 



— ;;2 ' 



< 



Pn Pn-1 

2Ci 



Pn-1 

in view of Ci > 2,p„ > p.„_i. Since |/3Q„„_i_i - i^«„_i-i| < Q^J_i, we have the same 
bound for \^nQu„-i~i - -P«„_i-i|, whence from (18.141 ). in view of |^^| < 1, 

r„-i < ^. (8.15) 

Pn-l 



By our construction, 






-Id n _ / Pn"+' ^ I Pn-1 - r-n-l 



whence 

\ '-' rn / \ • • • ' n-1 

Taking the traces yields 

Using r„ = Pn + Pn^ , we obtain 

Pn+l + Pnli = PTHPu-I + Pn-l ' ^n-l) + P;""+V„_i 

Therefore, 

1 _ P"+l _ ^«-l ,1 1 ^n-l 



rP-n+'^ n n n n°-n+l 2 2a„+i 

P«. Pn-1 Pn-1 Pn+lPn Pn-1 Pn-1 Pn Pn-1 

whence 

1 _ P?^+l ^ __}^ ^n-1 ^ C'o 



^a„+i ^ _ ^2 .2a„+i _ — ^2 



Pn Pn-1 P;-l pr"+>„_i Pn-l 

(in view of r„_i/p^""+i < 1 and Cq > 2), and 

1 _ Pn+l ^ ''"n-l I 1 



a„+i _ I a„+i 

Pn Pn-l Pn-l Pn+lpn Pn-l 

4Ci 



Pn-l Pn-l 

4Ci + 1 Co 



Pn-l Pn-l 



D 
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We also require the following lower estimate on the growth of the sequence (p^). 



Lemma 8.6. Ifn>l and K > 1 are integers such that the matrix i?n-i — K ■ I is 
non-negative, then pn+i > Kpn- In particular we have pn+i > 2pnfor all n > 3. 



Proof. Given a pair of matrices A and B we will use the notation A > B to mean that 
the difference A — B is a non-negative matrix. If A> B then obviously also tr A > tr B, 
and AC > BC and CA > CB for any non-negative matrix C. Note in particular that 
^0) ^1 > I^ and hence if C is any product of powers of Aq and Ai then C > /. It follows 
that Bn+i > Bn for all n > 0. 

lin.K > 1 and 5„_i > K ■ I, then B^+i = B^"+^Bn-i > KB^"+^ > KB^ and 
therefore r„_,_i > Krn. It follows that 

Pn+l = - [Tn+l + V(r„+i)2 - 4J > y f r„ + ^ (r„)2 - — J > /sfp^ 

as required. Since ai > 2, we may estimate 

B, = Bt'Bo = {B^,'-'B_^r^Bo = {A^,^~' A,)-' Ao 
> A,A,Ao =(l fj> 2/, 



and since Bn-i > -B2 for all n > 3 it follows that pn+i > 2pn for all n > 3 as claimed. 

D 



We may now give the proofs of Proposition l8.2l and Corollary 18.41 



Proof of Proposition \8. 2\ Let N > uq and define 
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Using (|8.6I) together with the second clause of Leninia l8.6l 



log a^ — log ajsf I 



J2 (-!)>« log ( 



Pn"+Vn-1 



log 



n=Af+l 

oo 

< 5^ g„ 

n=Af+l 

oo 

< Y. 9" 

n=Af+l 

oo 

n=Ar+l /^"-l 



Pn+1 



Pn+l 



Pn Pn-1 



Pn+l 



Pn Pn-1 



qn 



< 



oo ^ 



n=Af+l 
oo 

<^Co Y. 

n=N+l 

_ 2LCo 

Pn 



Pn-1 



2^+i-"Pjv 



as required. D 

Proof of Corollary \8A\ For each j > 1 we have 

Qno+l+j ^ ^no+l+jQno+j ' Qno+j-1 — -'^Qno+ji 

and it follows that qn^+i+j < -ft'-'gno+i for ^U J > 0. On the other hand since -B„(,-i — KI 
is non-negative, Bn„+j-2 is non-negative for every j > 1, and using Lemma 18.61 we 
deduce that Pno+j > -^Pno+j-i for all such j. We therefore have Qno+i+j < K^Qno+i < 
K^Lpno < Lpno+j for all j > and we may apply Proposition |8.2[ D 

Let us show that the hypothesis g^ = 0(p„_i) is valid for "typical" 7 in a suitable 
sense: 

Lemma 8.7. 7^7 is not Liouville, then qn < Pn-ifor all sufficiently large n. 

Proof. Since 7 is not Liouville, there exists 5 > such that 



Pn 

7 

qn 



> 



1 



« 



,5+1 • 



Since 



7 



Pn 

qn 



< 



QnQn+l 

we have g„+i < q^. Thus, it suffices to show that g^ < p„ for n large enough. By 
Lemma [831 pn ~ Pn-iPn-2, whence logp„ ~ a„logp„_i + logp„_2. Consequently, 
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logpn > const • g„. (Since g„ = a„g„_i + qn-2-) Now the claim follows from the fact that 
the Qn grow at least exponentially fast, whence log g„ ^ g„. D 

In fact the upper bound g„ = 0(p„_i) holds for "most" Liouville numbers as well. 
Effectively, if this inequality fails, this means that a„ > A'^"-^ infinitely often for some 
constant A > I, which is an exceptionally strong condition. 

Remark 8.8. It is natural to ask whether a formula like (18.21) - with traces instead of 
spectral radii - holds in a more general case of irrational 7 (instead of (12.21) . where the 
multipliers are irrational). The answer is yes - provided, for example, the condition g„ = 
O(pn-i) holds. Indeed, this condition implies 

1 \ Qn + l 

1 + ^r) -)-0 n-)-oo, 

PnJ 

whence we can replace the spectral radii with the corresponding traces so as to obtain 

«,= lim V =n ^^^^^ • (8.16) 

Note that if the a„ grow extremely fast (for instance, if a„ = g„_i)» then (18.161 ) is false; one 
can show that if it were true, then t^^ (7) would be an interval, contradicting Theorem l2.3l 
(iii). 

For 7 = ^~2 the formula (18.161) is exactly (18.21 ). in view of the recurrence rela- 
tion (|8.1I) . Indeed, we have g„ = F„ and 

n=0 ^ n=0 ^ ' 

Remark 8.9. Despite having such a fast convergent infinite product for a^, we still cannot 
use it to claim that a^ is irrational if 7 is irrational. Such a result would show that the 
family (11.11) does not contain a counterexample to the rational finiteness conjecture (see 
ifTTl for more detail). 

Remark 8.10. Another natural question is whether there exists a recurrence relation - or 
rather a sequence of such relations - for the r„ in the case of a general irrational 7. It can 
be shown that if a„ and a„+i are fixed, then there will be the same recurrence relation for 
r„+i, irrespective of the rest of a^. However, even in the simple case a„ = a„+i = 2, for 
instance, we have the relatively unstraightforward identity 

2 

n- — ^2^ "^n TnTn-2 

in+1 — 'n'n-1 'n-l- 

Tn-1 rn-l 

And for larger a„ and a„+i, it becomes messier, though the two most significant terms 
are always r^""+^r„_i — r^"+i/r„_i, provided the a„ do not grow too fast. The authors are 
grateful to Kevin Hare for helping them with these computations. 
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The following examples yield new explicit parameters a such that the system {Aq, a Ax } 
does not possess the finiteness property: 

Example 8.11. Put 7 = \fh — 2. It is algebraic and therefore, not Liouville. Here a = 

0.4596704785 . . . 

Example 8.12. Put7 = ^-1 = [3, 1,5, 1, 1,4, 1, 1,8, 1, . . .].Herea = 0.5587336687. . . 

Example 8.13. As is well known, e — 2 = [1,2, 1, 1, 4, 1, 1, 6, 1, 1,8,...], which implies 
that e is not Liouville. Put 

7 = ^^ = 0.4180232931... = [2,2,1,1,4,1,1,6,1,1,8,...]. 
e — 1 

Here a = 0.7904851693 . . . 
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